REPORT  DOCUMENTATION  PAGE 


Form  Approved  OMB  NO.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions, 
searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments 
regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggesstions  for  reducing  this  burden,  to  Washington 

Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington  VA,  22202-4302. 

Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  oenalty  for  failing  to  comply  with  a  collection  of 
information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY) 

2.  REPORT  TYPE 

3.  DATES  COVERED  (From  -  To) 

New  Reprint 

- 

4.  TITLE  AND  SUBTITLE 

Origin  of  second-harmonic  generation  enhancement  in  optical 
split-ring  resonators 


5a.  CONTRACT  NUMBER 
W91  INF-09-1-0539 


5b.  GRANT  NUMBER 


6.  AUTHORS 

Ekaterina  Poutrina,  Cristian  Ciraci,  Michael  Scalora,  David  R.  Smith 


5c.  PROGRAM  ELEMENT  NUMBER 

611103 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAMES  AND  ADDRESSES 

Duke  University 
2200  West  Main  Street 
Suite  710 

Durham,  NC  27705  -4010 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND 
ADDRESS(ES) 

U.S.  Army  Research  Office 
P.O.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 
ARO 


11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

561 54-PH-MUR.  85 


12.  DISTRIBUTION  AVAILIBILITY  STATEMENT 
Approved  for  public  release;  distribution  is  unlimited. 


13.  SUPPLEMENTARY  NOTES 

The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and  should  not  contrued  as  an  official  Department 
of  the  Army  position,  policy  or  decision,  unless  so  designated  by  other  documentation. 


14.  ABSTRACT 

We  present  a  study  of  the  second-order  nonlinear  optical  properties  of  metal-based  metamaterials.  A 
hydrodynamic  model  for  electronic  response  is  used,  in  which  nonlinear  surface  contributions  are  expressed  in 
terms  of  the  bulk  polarization.  The  model  is  in  good  agreement  with  published  experimental  results,  and  clarifies 
the  mechanisms  contributing  to  the  nonlinear  response.  In  particular,  we  show  that  the  reported  enhancement 
of  the  second  harmonic  in  split-ring  resonator  based  media  is  driven  by  the  electric  rather  than  the  magnetic 


15.  SUBJECT  TERMS 
Plasmonics 


16.  SECURITY  CLASSIFICATION  OF: 

a.  REPORT 

UU 

b.  ABSTRACT 

UU 

c.  THIS  PAGE 

UU 

17.  LIMITATION  OF 
ABSTRACT 

UU 


15.  NUMBER  19a.  NAME  OF  RESPONSIBLE  PERSON 

OF  PAGES  David  Smith _ 

19b.  TELEPHONE  NUMBER 

919-660-8258 


Standard  Form  298  (Rev  8/98) 
Prescribed  by  ANSI  Std.  Z39. 1 8 


Report  Title 

Origin  of  second-harmonic  generation  enhancement  in  optical  split-ring  resonators 

ABSTRACT 

We  present  a  study  of  the  second-order  nonlinear  optical  properties  of  metal-based  metamaterials.  A 
hydrodynamic  model  for  electronic  response  is  used,  in  which  nonlinear  surface  contributions  are  expressed  in 
terms  of  the  bulk  polarization.  The  model  is  in  good  agreement  with  published  experimental  results,  and  clarifies 
the  mechanisms  contributing  to  the  nonlinear  response.  In  particular,  we  show  that  the  reported  enhancement 
of  the  second  harmonic  in  split-ring  resonator  based  media  is  driven  by  the  electric  rather  than  the  magnetic 
properties  of  the  structure. 


REPORT  DOCUMENTATION  PAGE  (SF298) 
(Continuation  Sheet) 


Continuation  for  Block  13 


ARO  Report  Number  56154.85-PH-MUR 
Origin  of  second-harmonic  generation  enhancem 


Block  13:  Supplementary  Note 

©  2012  .  Published  in  Physical  Review  B,  Vol.  Ed.  0  85,  (20)  (2012),  ( (20).  DoD  Components  reserve  a  royalty-free, 
nonexclusive  and  irrevocable  right  to  reproduce,  publish,  or  otherwise  use  the  work  for  Federal  purposes,  and  to  authroize 
others  to  do  so  (DODGARS  §32.36).  The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and 
should  not  be  construed  as  an  official  Department  of  the  Army  position,  policy  or  decision,  unless  so  designated  by  other 
documentation. 


Approved  for  public  release;  distribution  is  unlimited. 


RAPID  COMMUNICATIONS 


PHYSICAL  REVIEW  B  85,  201403(R)  (2012) 

Origin  of  second-harmonic  generation  enhancement  in  optical  split-ring  resonators 

Cristian  Ciraci,1,  Ekaterina  Poutrina,1  Michael  Scalora,2  and  David  R.  Smith1 

lCenterfor  Metamaterials  and  Integrated  Plasmonics,  Duke  University,  Durham,  North  Carolina  27708,  USA 
2C.M.  Bowden  Research  Facility,  U.S.  Army,  Research,  Development  and  Engineering  Command  (RDECOM), 

Redstone  Arsenal,  Alabama  35803,  USA 

(Received  18  January  2012;  revised  manuscript  received  2  April  2012;  published  15  May  2012) 

We  present  a  study  of  the  second-order  nonlinear  optical  properties  of  metal-based  metamaterials.  A 
hydrodynamic  model  for  electronic  response  is  used,  in  which  nonlinear  surface  contributions  are  expressed  in 
terms  of  the  bulk  polarization.  The  model  is  in  good  agreement  with  published  experimental  results,  and  clarifies 
the  mechanisms  contributing  to  the  nonlinear  response.  In  particular,  we  show  that  the  reported  enhancement 
of  the  second  harmonic  in  split-ring  resonator  based  media  is  driven  by  the  electric  rather  than  the  magnetic 
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Metamaterials  (MMs)  are  artificially  structured  media 
whose  collective  electromagnetic  properties  derive  from  the 
geometry  of  sub  wavelength  inclusions.  To  date,  the  most 
common  MM  designs  have  made  use  of  inclusions  formed  by 
conducting  materials  that  function  as  subwavelength  electrical 
circuits.  These  conductor  based  MMs  have  been  proven  to  be 
adept  at  mimicking  a  wide  variety  of  linear  electromagnetic 
responses,  providing  a  unique  venue  to  explore  otherwise 
inaccessible  concepts.1  In  the  context  of  nonlinear  response, 
however,  artificial  materials  may  offer  even  greater  oppor¬ 
tunities,  due  to  the  inherently  inhomogeneous  local  field 
distribution  that  exists  within  and  around  MM  inclusions.  By 
carefully  structuring  the  inclusion  geometry,  extremely  large 
field  enhancement  regions  can  be  produced  that  can  dominate 
and  enhance  the  effective  nonlinear  response  of  the  composite. 

The  enhancement  of  nonlinear  processes  by  MMs  has 
been  demonstrated  at  radio  and  microwave  frequencies,  using 
packaged  components,  such  as  varactor  diodes,  to  introduce 
nonlinearity  into  the  gaps  of  metal  MM  inclusions.2  However, 
to  achieve  nonlinear  optical  materials  at  higher  wavelengths, 
a  simple  scaling  of  these  prototype  structures  to  higher 
frequencies  (e.g.,  beyond  a  few  terahertz)  will  not  suffice. 
First,  the  response  of  most  metals  changes  from  conductorlike 
to  dielectriclike  at  frequencies  above  a  few  terahertz,  with 
absorption  increasing  significantly  as  the  fields  are  able  to  pen¬ 
etrate  further  into  the  metal.  Second,  packaged  semiconductor 
components  are  not  readily  available  at  frequencies  above 
100  GHz. 

While  metals  and  conductors  may  possess  undesirable 
properties  at  optical  wavelengths,  such  as  increased  absorp¬ 
tion,  they  also  possess  unique  and  potentially  advantageous 
properties.  In  addition  to  large  field  enhancements,  metal 
nanostructures  also  support  intrinsic  nonlinearities  that  relate 
to  the  dynamics  of  free  and  bound  charge  carriers.  As  a  result, 
metals  possess  some  of  the  largest  nonlinear  susceptibilities 
known.  Examples  include  the  large  x values  of  gold  or  silver, 
for  example,  suggesting  that  metals  can  serve  both  to  form  the 
linear  MM  response  by  tailoring  the  structure,  while  serving 
as  the  source  of  nonlinearity  for  nonlinear  optical  MMs. 

The  second-order  nonlinearity  in  metals  arises  from  both 
volume  and  surface  contributions.  Nonlinear  surface  contri¬ 
butions  are  strictly  related  to  the  response  of  the  electrons 


within  the  Thomas-Fermi  screening  length  (A.tf^0.1  nm 
for  gold)  from  the  metal  boundaries.  In  this  subnanometer 
realm,  electron-electron  interactions  become  appreciable  and 
nonlocal  effects  must  be  taken  into  account.  Moreover,  since 
metals  are  centrosymmetric,  they  do  not  possess  an  inherent 
X(2)  nonlinearity.  However,  the  surface  of  a  metal  can  break 
spatial  symmetry  and  provide  a  mechanism  for  an  effective 
X(2)  nonlinearity.  This  homogenized  x(2)  nonlinear  response 
becomes  highly  dependent  on  the  metal  geometry,  making  it 
inherently  a  MM  construct. 

A  steady  stream  of  works  concerning  harmonic  gen¬ 
eration  from  metallic  nanostructures  has  been  published 
recently.3-7  Second-harmonic  generation  (SHG)  has  also  been 
observed  experimentally  from  a  variety  of  metal  nanoparticle 
systems,8-18  and  specifically  from  optical  split-ring  resonators 
(SRRs),  which  at  first  appeared  to  exhibit  an  anomalously 
high  conversion  efficiency.12,19,20  The  magnetic  resonances 
associate  with  SRRs  have  raised  speculation  that  the  enhanced 
SHG  results  from  a  strong  nonlinear  magnetic  response 
associated  with  the  Lorenz  force.19  Yet,  to  the  best  of  our 
knowledge,  a  convincing  explanation  of  the  nature  of  this 
conversion  enhancement  remains  lacking  in  the  literature. 

In  this  Rapid  Communication,  we  show  that  the  basic 
characteristics  of  nonlinear  optical  SRRs  may  be  explained 
solely  by  the  electric  properties  of  the  structure  rather  than 
its  magnetic  response,  in  contrast  to  the  conclusions  drawn 
from  previous  works.12,19  The  nonlinear  optical  response 
of  the  charge  carriers  in  the  metal  here  is  described  by  a 
hydrodynamic  model,  which  includes  the  effects  of  pressure 
associated  with  the  electron  gas.  To  facilitate  the  numerical 
models  and  remove  ambiguities  associated  with  additional 
boundary  conditions,  we  develop  an  expression  for  the 
effective  second-harmonic  surface  currents  in  terms  of  the 
polarization  vector  in  the  bulk  regions.  This  expression  allows 
us  to  easily  study  the  SHG  process  from  metal  particles  of 
arbitrary  shape. 

In  the  context  of  the  hydrodynamic  model,  the  electron 
fluid  density  n(r,t )  and  the  electron  velocity  field  v(r,0  satisfy 
Euler’s  equation, 

3v  e  B2 

- 1-  (v  •  V)  v  +  yy  =  —  (E  +  v  x  B) - Wn,  (1) 

3 1  m*  n 
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along  with  the  continuity  equation,  V  •  J  =  —eh,  with  J  = 
en\.  In  Eq.  (1),  m*  is  the  effective  electron  mass,  and  y 
is  the  electron  collision  rate.  The  last  term  in  the  equation 
is  due  to  the  electron  gas  pressure,  here  described  within 
the  Thomas-Fermi  model,  with  /3  proportional  to  the  Fermi 
velocity  vF.  Combining  the  continuity  equation  with  Eq.  (1), 
and  expanding  all  fields  in  a  perturbative  manner,  one  finds 
that  the  free  electron  polarization  P  =  J  satisfies  the  following 
set  of  inhomogeneous  equations:7 

2 

/S2V  (V  •  PO  +  (ft)2  +  icoyWx  =  (2a) 

m* 

£2V  (V  •  P2)  +  (4 «2  +  2/ )P2  =  -^E2  +  SNL,  (2b) 

m* 

where  the  nonlinear  source  Snl  is  given  by 

e  icoe 

SNL=  — Ei  (V  -P0+  —Pi  X  Bi 
m*  m* 

co2 

- [(V  •  P0  Pi  +  (Pi  •  V)Pi] .  (3) 

n0e 

The  subscripts  refer  to  fundamental  and  second-harmonic 
fields,  respectively.  In  deriving  Eqs.  (1)  a  harmonic  time  de¬ 
pendence  ( e~l0)t )  has  been  assumed.  Together  with  Maxwell’s 
equations,  these  equations  describe  the  fundamental  and 
second-harmonic  polarization  vectors,  respectively,  and  hold 
under  the  assumption  that  the  fundamental  field  remains 
undepleted.  The  nonlinear  source  given  by  Eq.  (3)  groups 
surface  and  bulk  second-harmonic  contributions.  Specifically 
we  have  the  nonlinear  Coulomb  term  (referred  to  as  a 
quadrupolelike  term  by  virtue  of  its  form)  proportional  to 
Ei  (V  •  Pi),  the  magnetic  Lorentz  force  contribution  Pi  x  Bi, 
and  convective  terms  (V  •  Pi)Pi  and  (Pi  •  V)Pi. 

The  effect  of  the  electron  gas  pressure  on  the  polarization 
vectors  Pi  and  P2  in  Eqs.  (2)  is  a  linear,  nonlocal  contribution 
of  the  form  f$2S/  (V  •  P).  This  term  has  been  predicted  to  be 
responsible  for  unusual,  resonantlike  phenomena  above  the 
plasma  frequency.21,22  In  general,  this  term  becomes  important 
in  a  region  of  order  of  Af  in  the  vicinity  of  the  metal  surface, 
where  electron-electron  interactions  dominate  the  nonlinear 
surface  sources  given  by  Eq.  (3).  However,  the  linear  nonlocal 
term  may  be  safely  neglected  in  Eq.  (2b),  since  it  does  not 
affect  the  amount  of  generated  harmonic.23 

The  presence  of  spatial  derivatives  (nonlocality)  in  the 
description  of  the  fundamental  polarization  vector  of  Eq.  (2a) 
requires  the  specifications  of  additional  boundary  conditions  to 
solve  the  electromagnetic  boundary  value  problem.24,25  From 
the  continuity  equation  and  Gauss’  theorem  one  obtains  that 
n  •  P  =  0  at  the  boundary.22,26  A  consequence  of  the  linear 
pressure  term  is  the  nonzero  extension  of  the  induced  surface 
charge  inside  the  metal.  That  is,  the  induced  electron  charge 
density  ^  V  •  Pi  is  zero  in  the  bulk  region,  and  rapidly  changes 
near  the  metal  surface,  where  it  reaches  its  maximum  value. 
This  behavior  reflects  the  variation  experienced  by  the  normal 
component  of  Pi ,  which  continuously  goes  to  zero  at  the  metal 
boundary. 

The  hydrodynamic  model  is  a  simplistic  model  of  electrons 
that  nevertheless  gives  a  fairly  accurate  description  of  linear 
and  nonlinear  processes  occurring  at  the  surface  of  metal¬ 


lic  structures.  However,  the  simultaneous  manifestation  of 
microscopic  and  macroscopic  scales  in  the  problem  makes 
the  resolution  of  the  system  of  (nonlocal  and  nonlinear)  equa¬ 
tions  quite  complex,  and  ordinarily  necessitates  considerable 
computational  resources  even  for  particles  whose  dimensions 
are  a  few  tens  of  nanometers.  Deep  inside  the  metal,  in 
the  bulk  region,  the  electron  pressure  may  be  neglected, 
since  it  leads  to  corrections  of  order  (A/tfA)2  l.23  On  the 
other  hand,  second-harmonic  conversion  efficiency  strongly 
depends  on  the  behavior  of  Pi  at  the  metal  surface  where 
the  impact  of  the  pressure  term  may  be  critical.  Moreover, 
if  the  electron  pressure  in  Eq.  (2a)  is  neglected  at  the  surface, 
the  theory  becomes  inherently  ambiguous,  as  pointed  out  by 
Sipe  et  al.23 

To  understand  where  the  ambiguity  arises  consider  the 
nonlinear  Coulomb  term  in  Eq.  (3),  which  is  proportional 
to  (V  •  Pi)Ei.  In  the  free  electron  limit  (/3  =  0)  the  normal 
component  of  the  electric  field  En  =  n  Ei  changes  discon¬ 
tinuous^  across  the  metal  boundaries.  This  discontinuity  is 
attributed  to  charge  accumulation  occurring  at  the  surface, 
which  assumes  the  form  of  a  Dirac  delta  function,  namely, 
V  •  Pi  =  Pn8o(n  •  r),  where  Pn  is  the  bulk  polarization  com¬ 
ponent  normal  to  the  surface.  It  follows  that  the  quantity 
(V  •  Pi )En  =  Pn8oEn  is  not  well  defined  and  cannot  be 
integrated.  The  same  analysis  applies  to  the  convective  terms. 
In  early  works  this  problem  was  circumvented  by  introducing 
phenomenological  coefficients  of  order  unity,27  or  through  an 
effective  plasma  frequency,23  that  incorporates  the  details  of 
the  charge  distribution  near  the  surface,  an  effect  that  was 
neglected  in  Ref.  6.  However,  an  analysis  of  the  effect  of 
pressure  on  the  amount  of  converted  energy  shows  that  the 
total  SHG  converges  to  an  asymptotic  value  as  f$  tends  toward 
zero.  We  verified  that  for  typical  values  of  f$  (of  order  of 
106  m/s),  the  actual  converted  energy  differs  by  not  more  than 
3%  from  its  asymptotic  value.  Exploring  the  limit  for  — >  0 
seems  then  a  very  good  way  to  get  an  approximate  solution 
without  having  to  solve  the  complex  nonlocal  equations,  as 
discussed  below. 

In  this  limit  the  surface  layer  is  so  small  compared  to  the 
size  of  the  nanoparticle  that  only  the  derivative  along  the 
direction  §  =  r  •  n  normal  to  the  metal  boundary  matters.  Let 
us  choose  the  metal  boundary  such  that  for  §  <0  the  electron 
pressure  may  be  neglected  (ft  =  0);  the  region  0  <  §  < 
where  /  ~  A.Tf,  represents  the  surface  layer  where  the  pressure 
is  important.  In  this  region,  the  parallel  component  of  the 
vector  Pi  is  nearly  constant,  pf  (§)  —  Pf  (0_),  while  its  normal 
component  P^~  may  be  written  as 

pt  (O  =  A(o_)°r(0>  (4) 


where  cr(i-)  is  an  unknown,  rapidly  varying  function  that 
goes  from  cr(0)  =  1  to  cr(Z)  =  0  (as  imposed  by  boundary 
conditions).  It  is  straightforward  to  show  that,  starting  with 
Eq.  (2b),  the  nonlinear  surface  polarization  is  given  by 


P 


s 

2 


1 

2noe 


(V  •  POP!  + 


CO 

2co  +  iy 


(Pl-V)P! 


(5) 


where  the  second-harmonic  pressure  term  /?2V(V  •  P2)  has 
been  neglected,  as  already  pointed  out.  In  writing  Eq.  (5)  we 
have  used  Eq.  (2a)  to  express  the  electric  field  Ei  as  a  function 
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of  Pi.  The  crucial  quantity  here  is  the  effective  nonlinear 
surface  current  density,  defined  as  Knl  =  2 ico  Pf  (r)dr, 

where  the  integral  is  performed  across  the  surface  layer.  Using 
Eq.  (5)  we  obtain 


KNl  = 


ICO 


n0e 


/v  i  II  1  3 co  ~\~  iy  19 
2  2  co  +  iy 


(6) 


where  cr'd%  =  —  1  and  crcr'd §  =  — 1/2  have  been  used, 
with  the  prime  representing  the  derivative  with  respect  to 
§,  and  the  unit  vector  t  pointing  in  the  direction  n  x  P25. 
Equation  (6)  states  that  surface  contributions  to  the  second- 
harmonic  polarization  may  be  approximated  by  an  effective 
nonlinear  current  sheet  at  the  surface  of  the  nanoparticle. 
These  currents  are  related  to  the  polarization  values  in  the  bulk 
region  and  do  not  require  the  resolution  of  nonlocal  equations. 
That  is,  the  volume  sources  are  calculated  by  assuming  ft  =  0 
while  the  surface  nonlinear  currents  are  given  by  our  Eq.  (6). 
In  the  limit  of  validity,  this  approach  provides  a  description 
of  the  SHG  process  that  may  be  easily  implemented  in  full 
three-dimensional  (3D)  simulations.  We  tested  this  approach 
in  the  case  of  SHG  from  a  silver  slab  and  compared  our  results 
with  experimental  data.28  The  results  obtained  using  Eq.  (6) 
agreed  very  well  with  both  the  solution  given  by  the  complete 
Eqs.  (1)  and  the  experimental  data. 

Three-dimensional  metal  nanoparticle  systems  are  of  con¬ 
siderable  interest  for  nonlinear  media,  because  both  the  field 
enhancement  regions  and  the  surface  morphology  strongly 
impact  SHG.  In  an  experiment,  planar  arrays  of  varying  MM 
inclusions  were  shown  to  produce  second-harmonic  light, 
with  efficiencies  that  varied  according  to  the  shape  of  the 
nanostructured  inclusions.12  To  investigate  the  mechanism  of 
SHG  from  metal  nanoparticles,  we  choose  the  same  types  of 
nanoparticles  studied  in  Ref.  12.  Our  results  are  summarized 
in  Fig.  1  for  the  different  kinds  of  nanoparticles.  To  ease  the 
computational  burden  and  to  simplify  the  system  studied,  we 
assumed  the  nanoparticles  to  be  surrounded  by  air  rather  than 
including  the  substrate  material  used  in  the  experiments.  We 
do  not  expect  significant  changes  in  the  observed  mechanism 
of  SHG  due  to  the  exclusion  of  the  substrate  material.  The 
structure  studied  extends  periodically  in  the  x  and  y  directions, 
so  that  only  a  single  unit  cell  is  needed  in  the  computational 
space.  To  avoid  possible  numerical  artifacts  due  to  the 
field  localization  near  metal  comers,  we  considered  rounded 
corner  geometries  with  a  radius  of  curvature  of  5  nm.  The 
geometrical  parameters  were  chosen  so  that  the  nanoparticles 
would  display  a  resonance  around  Aff  =  1.5  /xm,  to  which 
the  fundamental  field  is  tuned.  The  conversion  efficiencies  r) 
assume  an  average  pump  intensity  of  ~55  MW/cm2  with  the 
electric  pumping  field  polarized  along  the  x  direction. 

We  find  qualitatively  good  agreement  with  the  experimental 
data  of  Ref.  12,  where  second-  and  third-harmonic  generation 
were  experimentally  investigated  for  rectangular,  T-  and  U- 
shaped  gold  nanoparticles,  respectively.  We  find  the  relative 
efficiencies  for  normal  incidence  normalized  with  respect 
to  the  U-shaped  nanoparticle  efficiency  to  be  about  1.7% 
for  T-shaped  nanoparticles  (compared  with  2.1%  for  the 
experimental  data),  while  the  rectangular  nanoparticle  showed 
a  conversion  efficiency,  as  one  might  expect  due  to  the 
absence  of  symmetry  breaking. 
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FIG.  1.  (Color  online)  Second-harmonic  conversion  efficiency 
for  different  3D  gold  nanoparticles.  All  the  particles  are  20  nm 
thick.  The  linear  transmittance  (solid  line)  and  reflectance  (dashed 
line)  at  normal  incidence  are  shown  in  the  second  column  for  the 
variety  of  nanoparticles.  The  vertical  lines  indicate  the  fundamental 
(dotted  red  line),  and  the  second-harmonic  (blue  dashed-dotted 
line)  wavelength,  respectively.  In  the  third  column,  the  arrows 
denote  the  polarization  state.  The  fundamental  field  is  represented 
in  red  (light  gray),  and  the  second-harmonic  field  is  blue  (gray). 
The  relative  conversion  efficiencies  normalized  with  respect  to  the 
U-shaped  nanoparticle  efficiency  are  also  reported  (the  corresponding 
experiment  of  Ref.  12  are  shown  inside  brackets).  The  following 
values  for  the  parameters  have  been  used:  ra*  =  me,  n0  =  5.7  x  1022 
cm-3,  and  y  —  1.07  x  1014  s-1. 


As  in  Ref.  12  we  find  that  an  array  of  U-shaped  gold 
nanoparticles  can  enhance  the  SHG  efficiency  by  about  two 
orders  of  magnitude  with  respect  to  other  noncentrosymmetric 
nanoparticles,  such  as  T-shaped  particles,  for  example.  More¬ 
over,  the  polarization  state  of  the  generated  field  is  rotated  by 
90°  for  a  pumping  electric  field  polarized  along  the  x  direction. 
This  polarization  flip  is  not  surprising,  because  symmetry 
breaking  occurs  along  the  y  direction.  However,  the  nature 
of  the  conversion  enhancement  remains  an  unsettled  point 
in  the  literature.  In  Ref.  19  the  authors  speculated  that  this 
effect  arose  as  a  result  of  a  strong  nonlinear  magnetic  response 
associated  with  the  Lorenz  force.  However,  the  Lorentz  force, 
proportional  to  Pi  x  Bi,  is  but  one  of  several  influential 
contributions,  Eq.  (3).  Within  the  context  of  our  model,  by 
simply  switching  off  the  Lorentz  force  contribution  to  the 
generated  field  we  find  that  the  conversion  efficiency  is  only 
marginally  affected.  The  effect  is  shown  in  Fig.  2(a),  where 
the  second-harmonic  conversion  efficiency  is  calculated  as  a 
function  of  the  pumping  wavelength  Aff,  with  and  without 
the  Lorentz  force  term.  Therefore,  the  magnetic  response  of 
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FIG.  2.  (Color  online)  (a)  Comparison  of  the  second-harmonic 
conversion  efficiency  with  and  without  Lorentz  force  contributions 
for  the  U-shaped  nanoparticle  of  Fig.  1.  (b)  Sketch  of  the  fundamental 
and  second-harmonic  currents  flowing  in  the  rectangular  and  U- 
shaped  nanoparticles,  respectively.  For  the  symmetric  nanoparticle 
the  nonlinear  currents  cancel  out  and  no  second-harmonic  field  is 
produced.  Instead,  for  the  U-shaped  nanoparticle  the  currents  oscillate 
in  phase  and  radiate  to  the  far  field. 


nanoparticle.  If  the  Lorentz  contribution  is  neglected,  all  the 
electric  nonlinear  contributions  are  basically  proportional  to 
the  product  of  the  electric  field  E  and  the  induced  charge 
distribution  p  =  —  en  =  V  •  P,  as  may  easily  be  deduced  from 
Eq.  (2b).  In  Fig.  2(b)  two  parallel  situations  are  depicted: 
(i)  the  rectangular  nanoparticle  and  (ii)  the  U-shaped  nanopar¬ 
ticle.  The  arrows  indicate  the  electric  field  at  the  surface, 
while  the  color  map  represents  the  charge  distribution.  At 
the  fundamental  frequency  both  structures  behave  similarly. 
The  external  electric  field  of  the  electromagnetic  radiation 
produces  a  gradient  in  the  electrical  potential  that  causes 
charges  to  migrate  from  one  pole  to  the  other.  To  a  first 
approximation,  the  resulting  second-harmonic  polarization  is 
proportional  to  the  product  between  the  charge  distribution 
and  the  electric  field,  that  is,  P(2 co)  oc  p(co)E(co).  As  shown 
in  the  sketch  of  Fig.  2(b),  the  nonlinear  polarization  will 
result  in  a  second-harmonic  electric  field  that  cancels  out 
for  the  centrosymmetric  nanoparticle,  resulting  in  a  near-zero 
second-harmonic  efficiency.  On  the  other  hand,  for  the  U- 
shaped  nanoparticle  the  resulting  harmonic  fields  oscillate  in 
phase  and  the  produced  field  can  radiate  to  the  far  field.  These 
findings  lead  us  to  suggest  that  the  characteristic  emission  of 
nonlinear  optical  SRRs  is  driven  by  the  electric  properties  of 
the  structure  rather  than  its  magnetic  response,  as  originally 
conjectured. 

In  conclusion,  we  have  discussed  the  influence  of  the 
electron  gas  pressure  on  SHG  and  find  a  way  to  express 
nonlinear  surface  contributions  in  terms  of  the  polarization 
vector  in  bulk  regions.  Using  this  approach  we  have  analyzed 
the  SHG  in  several  types  of  3D  gold  nanoparticles  and 
demonstrated  that  the  origin  of  the  previously  experimentally 
reported  enhanced  SHG  from  SRR-based  nanostructures  is 
mostly  electrically  driven.  This  development  should  help 
simplify  the  investigation  of  the  SHG  process  in  full  3D 
metal  structures,  and  enable  the  investigation  of  arbitrary 
nanostructured  plasmonic  media,  offering  a  viable  design 
approach  to  integrated,  efficient,  nonlinear  optical  media. 


the  U-shaped  nanoparticle  does  not  seem  to  play  an  important 
role  in  the  description  of  the  SHG  process,  at  least  assuming 
the  hydrodynamic  model. 

The  origin  of  the  conversion  enhancement  may  be  ex¬ 
plained  by  studying  the  nonlinear  electrical  properties  of  the 
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